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, Abstract 

b ■ 

, In this paper we discuss the Mather problem for stationary Lagrangians, that is La- 

grangians L : R" x R" x f2 ^ R, where f2 is a compact metric space on which R" acts through 
an action which leaves L invariant. This setting allow us to generalize the standard Mather 
0^ ■ problem for quasi-periodic and almost-periodic Lagrangians. Our main result is the existence 

of stationary Mather measures invariant under the Euler-Lagrange flow which are supported 
in a graph. We also obtain several estimates for viscosity solutions of Hamilton- Jacobi equa- 
tions for the discounted cost infinite horizon problem. 
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1 Introduction 

Let M be a complete compact manifold, and L : TM ^ R a Lagrangian, fiberwise strictly 
convex and coercive. A probability measure on TM is called holonomic if 

'nT , Jtm 

0^ , . 

lO ■ for all ip e C^(Af). A central result in Aubry-Mather theory [Mn96j (see also |FS04| ). is the 

fact that any holonomic probability measure ^ on TM which minimizes the action Jrpj^j Ld^ is 
CO ' supported on a Lipschitz graph and is invariant under the Euler-Lagrange flow. Certain results 

I in Aubry-Mather theory have been extended for non-compact manifolds, see for instance |FM07j 

Q*^ ■ or [Mad06j . but as far as the authors know, there is in the literature no satisfactory construction 

of Mather measures for general non-compact manifolds. 

In this paper, rather than considering Lagrangians on the tangent bundle of compact manifolds, 
such as in the original paper of Mather [Mat91| . we consider Lagrangians defined on K" x R" x il, 
; ^ ' where f2 is a suitable compact metric space on which R" acts trough an action . The main result 

5^ \ of this paper is Theorem in which we establish the existence of stationary Mather measures 

invariant under the Euler-Lagrange flow. 

Stationary ergodic problems were considered in [LS03| in the context of homogenization of ran- 
dom stationary ergodic Hamilton- Jacobi equations. The authors (in particular DG) are thankful 
to several enlightening discussions with P. Souganidis on this issue. Generalized Mather measures 
for stationary ergodic problems were also considered in the homogenization setting in {GV07| . 
The stationary ergodic setting was consider in [?] where the construction of critical (or criti- 
cal approximate) viscosity solutions of Hamilton- Jacobi equations is carried out in detail for the 
onc-dimensional case. 
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A simple example (taken from |LS03j ) which illustrates the main difficulties in the stationary 
setting is the Lagrangian 

L = — cos(x + uji) — cos(\/2a; + uj2)- 

Consider uj E M^/Z^ = as a fixed parameter. It would be natural, as in Mather's problem, to 
look for probability measures /i on M" x M" which minimize the action 

Ld^l (1) 

under the holonomy constraint 

V ■ Dxtpdfi = 0, 



for all (p of class C^, bounded with bounded derivatives. This problem can be solved explicitly, 
and in fact we have the following two cases: if there exists a solution x to the overdetermined 
system 

X + LOi = 27rn, V2x + Ll>2 ^ 27rn, 

for some n € Z, the Mather measure on M x M is simply = dx{x)5Q{v); otherwise there does 
not exist a Mather measure since L > —1 for all and the infimum in ([1]) is easily shown to 

be -1. 

To overcome these issues, which are due to the lack of compactness of M", we will instead 
define stationary Mather measures as measures on {v,uj) € R" x fl, which minimize the action 
and satisfy a suitable holonomy condition. It turns out that if f2 is compact and the Lagrangian 
satisfies certain stationarity hypothesis this is the natural way to generalize Mather measures. 
Before proceeding, we must make precise our framework. 

Let 17 be a compact metric space, and let L = L{x, v,uj) : R" x R" x i7 R be a continuous 
Lagrangian, in the first two coordinates. The Lagrangian L is also required to be strictly 
convex and superlinear on the velocity v, and nonnegative. In our setting, this last condition can 
be achieved without changing of the nature the problem by adding a constant to L. We assume 
further that 

L{x + y, V, u) - L{x, V, < \y\ {C + CL{x, v, u)) . (2) 

We suppose that there exists an action r : x M" ft which is continuous, satisfies the 
semigroup property 

T^+yU; = T^TyLo and to(-) = Id. 

Since Vl is compact and the action is continuous, the action is uniformly transitive^ in the following 
sense: 

Ve > 0, 3M > 0,Vwi,cj2 e fi, 3z e M", such that \z\ < M, and d{T^uJi,LLi2) < £■ 

A first example of such an action is the following: we take 17 = T'*, the d-dimensional torus, 
let n < d and we will construct an action r : R" x T'' — > T"*. To start with, we identify the torus 
T"* with its universal covering K'', and consider a constant coefficient d x n matrix A. Assume 
that {Ax : x € R"} is dense in T". Then we define 

TxiV — LU + Ax. 

A second example is the following. We take to be the space of all sequences w = {ojk) on 
T^, endowed with the following metric: 

oc 
k=l 



^The authors are grateful to Albert Fathi that pointed out to us that uniform transitivity holds under the 
compactness assumption. 
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It is simple to verify that with this distance the space f2 is compact. A sequence A of real numbers 
is called irrational if for any N the vector (Ai, . . . , Ajv) is is linearly independent over the integers. 
Let A be an irrational sequence. Define the following action from M into H. by 

TxU> — UJ + xX. 

This action is also uniformly transitive. 

A function : R" x M" x — + R, is stationary if 

(p{x + y, V, oj) — (p{x, V, Ty{uj)), Vx, y £ M", a; g fi. 

We assume that the Lagrangian L is stationary. 
Denote 

C^(]R" X ri) = {(^ : M" X ^ M, stationary, in the first variable, continuous in 
OJ, and such that I?a;(^(0,Lj) is continuous in cj}, 

with an analogous definition for Cj(E" x R" x fi). 

If the action is given as in the first example by t^uj = o; + Ax, given ■(/; : T'^ ^ M, the function 
ip{x,Lo) = ip{u! + Ax) is stationary, and, furthermore, tp £ Cl if ip is . In the second example 
we can construct an example of a stationary function in the following way: let 4'k : T ^ M be a 
sequence of periodic functions uniformly bounded in k. Let 

fc l+|Afc| 

Furthermore, if ipk is and its derivatives are uniformly bounded in k, ip G Cl . 

To motivate the stationary Mather problem, let x{t) be a globally Lipschitz trajectory on M". 
Let Wo G is an arbitrary point. Consider ergodic averages to define an occupation measure /x 
on R" X il corresponding to x{-) in the following way 

lim ^ I 4i{x,x,uoQ)dt ^ lim \; [ (j){Q,x,TxU}Q)dt = [ (l){0,v,uj)dfj., 

where the limit is taken trough an appropriate sequence. Of course, the measure /i could depend 
on the point luq or the sequence through which the limit is taken. Nevertheless, such probabilities 
H, satisfy an integral constraint, the holonomy condition: 



V DMO,i^)d^i = 0, (3) 



"xSl 



for any stationary function ip e Cj(R" x fl). 

The stationary Mather problem can be formulated as follows: minimize 

L{0,v,uj)dfi{v,uj), 

'Xf2 

over all probability measures that satisfy the holonomy constraint ([3]). A minimizing measure for 
this problem is called a stationary Mather measure. A similar problem arises also in [GV07| for 
the homogenization of Hamilton-Jacobi equations. 
Let 7 : M" ^ M be a positive function such that 

lim = 0, and hm ^''^'^[^^ +00, (4) 

\v\^oo 7(v) |-u|^oo 7(f) 
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where the last hmit is uniform in a; e by compactness. We denote by C°(R" x il) the set of the 
continuous functions (f> with 

Uh = sup — — < oo, hm — > 0. 

We will need also to consider the discounted Mather problem, see |Gom08] for a discussion of 
related generalizations of Mather's problem. For that, let a be a positive number. Consider the 
operator A : Cl{W x n) ^ C°(R" x Q) given by 

ip A'"(p{ll>) — V ■ Dx(p{0, Lo) — a^{0, uj). 

The discounted stationary Mather problem consists in minimizing 

L{0, V, Lo)dp,{v, w) 

over all probability measures that satisfy the discounted holonomy constraint 

ip{uj)diM{v,uj) = —a I ip{Q^ULi)di'[ijj), (5) 



for all Lp e (R" xVt). A minimizing probability measure for this problem is called a discounted 
stationary Mather measure. The measure v is called the trace of /.t. If a = we call these measures 
stationary Mather measures. 

The main result of this paper is the construction of stationary Mather measures invariant 
under the Euler-Lagrange flow. Usually, this flow is defined in R" x R". However, since the 
stationary Mather measures are measures on R" x £7 we must now discuss the natural extension 
of the Euler-Lagrange flow to this space. 

Given a stationary vector fleld : R"xR"xl] ^ R"xR", let $ = ($i,$2) : RxR"xR"xl] ^ 
R" X R" be its flow. We define the flow : R x R" x O ^ R" x 17 induced by in R" x 17 as 

■^{t,V,Uj) = ($2(i, 0,-U,CJ),T<I,J(i^o,t;,c^)^)■ 
We denote by C^(R" x Vl) the set of bounded continuous functions 0(w,cj) in R" x Q, such that 
Dy(j){v,uj) is also continuous and bounded. A measure is invariant under the flow if, 



'xO JM"xQ 

for aU (f) G Cbi(R" x fl) and for aU t e R. 

Let fihe a measure in R" x fl and : R" x R" x il ^ R" x R" be a stationary vector field 
in R" X R". Then /x is invariant under the flow induced by W in R" x 17, if and only if, 

V(^(0, V, uj) ■ W{0, V, uj)dfi{v, iv) = 0, (6) 

'xO 

(where the gradient in the previous formula is taken both in x and v) for all (j) e (^^(R" x R" x 17). 
A proof for this classical fact for the case of vector fields on a manifold M can be found, for 
instance, in [BG08j . The proof in our setting follows exactly along the same lines and we will omit 
it. 

In this paper we will need to consider the discounted Lagrangian La = e^"*L(x,w). The 
corresponding Euler-Lagrange equation is 

^DyL{x,v,uj) ^ DxL{x,v,uj) + aDyL, (7) 
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for each uj Cz fl. For a = we obtain the usual Euler-Lagrange equations. We have a w-parametric 
Lagrangian vector field , that is given by: 




X'"" {x, v,uj) = V 



We say that a measure /i in K" x is invariant under the Euler-Lagrange flow if it is invariant 
under the flow vE"" induced by W^" in R" x 17. 

The outline of this paper is as follows: in section [3 we describe briefly the duality theory for 
the stationary Mather problem and its connections with viscosity solutions of Hamilton- Jacobi 
equations. The proofs of some the results, since they are standard, are outlined for completeness 
in appendix El In section [3] we make some formal computations in the spirit of [EGOl] . These 
computations suggest that for certain discounted stationary Mather measures one may be able to 
extend the regularity results in [EG01| . Holonomic discounted stationary Mather measures are 
constructed in section 3] Using these measures we obtain regularity results for viscosity solutions 
in section [51 These imply that the discounted stationary Mather measures are supported in 
a (partially) Lipschitz graph whose Lipschitz constant is independent of the discount factor a. 
Finally in the last section we construct stationary Mather measures invariant under the Euler- 
Lagrange flow. 



2 Duality and viscosity solutions 

The stationary Mather problem is an infinite dimensional linear programming problem. As usual 
in these problems (see jGomOSj . for instance), the duality theory plays an important role and will 
be developed in this section. 

Theorem 1. Let v he & probability measure on 51 and a > 0. Define 

Ha:= vni I L{0,v,uj)dii{v,uj), (8) 
Jr" X n 

where the infimum is taken over all probability measures on K" x O which satisfy the discounted 
holonomy condition Let 

Ti.{(p, X, uj) = sup {—A^{(p){x, Lo) — L{x, v, lo)) = -ff (x, D^^pix, w), uj) + aip{x, w), 

with H{x,p, u) = sup„gR„ {-p-v- L{x, v,uj)). 

Then, the infimum in ([8]) is achieved at some probability measure fi satisfying ([5|) and further- 
more 

Ha = — inf sup < — a / (pdv + Ti.(ip,0,Lu) > . (9) 

The proof of this Theorem is similar to analogous results in |Gom08j . for instance. For com- 
pleteness, however, we present the proof in the Appendix \X[ 

In this paper we will need to consider viscosity solutions to the equation 

H°(u,0,w) = i7(0,L»a.Ua(0,w),w)-t-aMa(0,w) =0. (10) 

As in the standard Mather problem, viscosity solutions yield important information concerning 
the value of the variational problem ([H]), and help characterize the support of the measure. 

Before we proceed, we make some remarks concerning the regularization by convolution of 
stationary functions. 
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Remark 1. To approximate a stationary function u : M" x il — > R by smooth stationary functions 

we arc going to use a convolution with a standard moUifier 77^ : M" R, that is. 7] compactly 
supported, ri^{x) = 7?7(f ), and /g„ r}{x)dx = 1. We define the convolution between u and ry*^ by 

u^{x,uj)= / u{x,TyUj)r]^{y)dy. 

Observe that, G C^. Moreover, we have 

du'^ f 

— {x,uj)-v = - J u{x,Ty(j)Dyr]%y) -v dy. 



We consider two different types of viscosity solutions for 0,a;) = A. Firstly recall the 
usual definition of viscosity solution: a function u : K" x 17 — !■ M, continuous in x (not necessarily 
C^) for each u £ il, is a. viscosity solution in x of T-L^u.x.lj) = A if for each loq e fi, any 
function t/' : R" — > R and any xq € M" such that u{x, loq) — i){x) has a strict local minimum (resp. 
maximum) at xq with u{xo,loo) — iI}{xq) = we have 

W(V',a;o, Wo) > A (resp. < A). 

For our purposes we need a modified version of viscosity solution: a stationary (not necessarily 
C^) function u : M" x J7 ^ R, continuous in f7, is a viscosity solution in uj of 7i('u, 0, w) = A if for 
any ip (R" x fi) and any point <= such that w(0, w) — i^(0, w) has a local minimum (resp. 
maximum) at with u(0, wq) — (/^(O, ujq) = we have 

H(93, 0,wo) > ^ (resp. < A). 

Proposition 2. Suppose that u : R" x f2 ^ R is a viscosity solution in x of H(w,0,a;) = A and 
assume furthermore that u is stationary and continuous in O. Then u is also a viscosity solution 
in w of W(u, 0, w) = A. 

Proof. Let ?i : R" x fi ^ R be a viscosity solution in x of TC{u,0,uj) = A. Consider an arbitrary 
fimction ip E (K" x fl) and a point o^o G such that u(0, w) — ip{0., ui) has a local minimum (resp. 
maximum) and m(0,wo) — </?(0,<x'o) = 0. Define ^{x) = (p{x,oJo)- We claim that u{x,(jJo) — ipix) 
has a local minimum (resp. maximum) in a;o = e R". In fact, 

u{x,U>o) - 1p{x) = u{x,OJo) - ip{x,Wo) = u{0,Ta;U)o) - 95(0, TxWo) 

> ^(0, Wo) - <^(0, Wo) = ^(0, Wo) - V'(O), (resp. < .) 
Then, because u is a viscosity solution in x we have 

H(V',0,wo) = H((p,0,wo) > A (resp. < A). 

□ 

Consider the infinite horizon optimal control problem 

Ua{x,u)) = inf / e~°'*L{x{t),x{t),uj)dt, (11) 

x(0)=xJq 

where the infimum is taken over all globally Lipschitz trajectories with initial condition x(0) = x. 
Then Ua ■ R" x O — > R satisfies the dynamic programing principle 

Ua{x,u})= inf 1 / e-"*L{x{t),x{t),u})dt + e-"'^Ua{x{T),u)y\ , (12) 

x{0)=x \Jq I 
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among all globally Lipschitz trajectories with initial condition a;(0) — x. It is standard, see 
[BCD97) . that the function Ua is a viscosity solution of 7i((/J, 0,1^) = in x. Furthermore, the 
optimal trajectories are solutions to the discounted Euler-Lagrange equations ([7]). Finally, for 
< < < r we have additionally that DxUa{x{t)) exists and 

xit) = -DpH{D,Uo,{x{t)),x{t)). 

The next proposition is also a well known result, see, for instance, |BCD97| for similar results: 

Proposition 3. For each oj fixed, let Ua{x,uj) be a viscosity solution (in x) of 

H"{u, a;, uj) — H{x, DxUa{x, uj), uj) + aUa{x, uj) = 0. (13) 

Then aua is uniformly bounded and Ua is uniformly Lipschitz in x, as a — > 0. 

Using standard techniques we can establish the following proposition, whose proof is presented 
in appendix [Bl 

Proposition 4. Let Ua '■ R" x ^ M be a solution of (fTS]) . Then Ua is a viscosity solution (in lo) 
of 7Y(i^,0, = 0, and Mq(0, w) is Lipschitz in w with Lipschitz constant (in uj) bounded by K/a, 
where K is independent of a, for all a > 0. 

Proposition 5. Let Ua be a viscosity solution in lo of (fT^l) Then 



inf sup < — a / ip{Q,u!)dv{uj) +'H°'{lp,Q,uj) \ — —a I Ua{0,Lu)di'{iL;). 
veci^gn [ Jn } Jn 

Proof. Consider a viscosity solution Ua of (fT^ . Then for any tp £ there exists a point uj^ of 
minimum for Uq,(0,w) — ip(0,uj). Consider ip'{x,uj) = ip{x,uj) + {ua — (p)(0, (jJj^). Then Ua{0,u!) — 
(/3'(0,a;) has a minimum equal to in uj^. 

Since Ua is a viscosity solution we have T-L°'{if' , 0, uj^^) > or equivalently 

n°'{^, 0, oj^) + a{u„ - (p)(0, oj^) > 0. 

Therefore 

-a (p{0,LLi)di^ + 'H°'{(p,0,uj^) + a{ua - (p){0,uj^) > -a (p{0,uj)di^, 
Jn Jn 

which implies 



and so 



sup —a j (y9(0, a;)(iz^ + ?i"((/5, 0, w) > — a j ip{0,uj)dv + a{ua — ^)(0,LUtp), 
wen 



sup —a / (p{0,uj)dh' + H°'{ip,0,Lu) > —a j Ua{0,Lu)di', 
uien 



which finally yields 

inf sup < — a / ip{0,uj)di^{Lj) +T-l"{ip,0,uj)> > —a / Ua{0,uj)di'. 
vec-i^^en I Jn } Jn 

In order to get the other inequality we use the functions = Ua Then 7i"(M^, 0, lo) < o(l) 
owing to the convexity of the Hamiltonian and the uniform Lipschitz estimates on Uq., we have 



inf snp<-a / ip{0,Lo)du{Lo) +n°'{ip,0,Lo) } < o{l) ~ a / u^O,Lo)diy. 
veci^en L Jn J Jn 

Then, the inequality desired is obtained by sending e to 0, and ends the proof. □ 
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Corollary 6. We have 

Ha = a I Ua{0,u})diy 



n 



where Ua is the unique viscosity solution of H{0, DxUa{0, + Q!Uq,(0, lo) = 0. 

Proof. In fact, if we apply Proposition [5] we have the formula 

inf sup < — a / ip{0,u!)dh'{uj) +Ti°'{(p,0,uj) > = —a / Wq(0, w)dz/(ti;). 
Remembering that Ha = — infi^gci sup^^gj-^ —a (p{0, uj)dv{uj) + T-l"{(p, 0, lo), we get 

Ha — Ua{0,Oj)dl'. 

Jn 



□ 



We state next, without proof, a partial converse to Proposition^ The proof is rather technical 
and, in this paper, its only a application is in Remark [51 

Proposition 7. Suppose that, 

(a) There exists S > such that, for all a; 7^ with < 5, Ta;(-) : ^ ft does not have fixed 
points. 

(b) For each luq G fi, there exists 5 > and a set 'Ssi'^o) 9 wq, such that, for all x 7^ with 
|a;| < S, and wi, (jJ2 G ^si^o), if Tx{uji) ~ L02 then loi — lo2- 

(c) The set 

Usiujo) = {t^{uj)\ uo e I],-(a;o), |a:| < 5/2} (14) 
is an open neighborhood of loq. 

If u : R" X fi ^ M is a viscosity solution in w of 7i (w, 0, w) = A then u is also a viscosity solution 
in X of 7^"(u,0,w) = A. 

Remark 2. Note that in some cases T-1°'{u,lo) = A does not admit viscosity solutions in w, 
as pointed out in [LS03j . In their example = T^, i = L{x,v.,uj) : R x R x ^ R is 
the Lagrangian given by L{x,v,uj) — + cos(ti;i + x) + cos(ti'2 + VSa^), with the associated 
Hamiltonian H{x,p,uj) = ^p^ — cos(a;i + x) — cos(ti>2 + V^x), and the action r : R x ^ T-^ is 
given by Tx{uji,uj2) = (wi +x,U!2 + V^x). 

In this case the viscosity solutions in x are unbounded. So, if there where a viscosity solution in 
to, then it would be a solution in x by Proposition [71 By compactness, any stationary continuous 
function is bounded, which would be a contradiction. 

3 Some formal computations 

In this section we adapt the formal computations in fEGOlj to motivate the regularity results in 
the following sections. Consider the periodic case of a Lagrangian L : T x R ^ R, given by 
L{x, v) = — V{x)^ and the associated Hamiltonian H{x,p) — ^p^ + V{x). The stationary case 
follows along the same lines, as we will see in later sections. 

Let M be a solution to the discounted Hamilton- Jacobi equation + V{x) + au = 0. Let 
fj,a be a discounted Mather measure with trace 9a and such that the projection of /^^ in the x 
coordinated is denoted by 9, that is, 

(p{x)d^a = I ^{x)d9. 

TxK Jt 
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Note that 9 in general does not agree with 6 a- In this section we assume that /x^ has the special 
property that 6^ = 0. Under this assumption fi^ is holonomic, that is 



/ Vipa:{x)dlJ,a = 0, 

JtxM. 



for all periodic function (f{x). 

We will first show that almost every {x,v) e T x M, we have v = —Ux{x). To see this we 
will argue by contradiction. In this case if t; ^ —Ux{x), there would exist a set of positive measure 
fia in which 

L{x, v) + vux > -H{ux, x). 
Since L{x,v) +vUx > —H{ux,x), integrating with respect to yields 

/ LdHa + / VUxdHa > a udna- 

This would yield 

/ LdiJ,a > ct ud9a, 

JtxR Jt 

which contradicts the optimality condition. 

Therefore the holonomy constraint can be written as 



/ {Ux^x 

Jm. 



aip)d9{x) = a I (fd9a{x). 



By differentiating twice the Hamilton- Jacobi equation we have Ux{uxx)x+Uxx+V' {^)+'^'^xx = 
0. Integrating with respect to yields 

+ ulx + V"{x) + aUxx)dfia = 0, 

Jm 

or, equivalently, 

/ Ux{uxx)xd9{x) + / ulx + auxxd9{x) = - / V"{x)d9{x). 

Since the trace of fia, 9a is equal to its projection ^(a;), then the measure fia is holonomic and so 

/ Ux{uxx)xd9{x) = 0. 
Jr 

Using —auxx < \u^x + 5^^ S^t, 

/ ulxd9{x) = - I V"{x)d9{x) - j auxxd9{x) 
Jr Jr Jr 

<- [ V"{x)d9{x)+ [ lulx + la'd9{x), 
Jr Jr ^ 

which yields a L'^{9) bound for Uxx' 

I ulxd9{x) < [ a^- 2V"{x)d9{x). 
Jr Jr 

In order to derive L°° estimates to Uxx we proceed as follows: first we multiply the second derivative 
of the Hamilton- Jacobi equation by a function ^''(uxx), 

/ '^xi'^xx^x 

'^'iuxx)d9{x) + / [ulx + aUxx + V"{x)\ ^'{uxx)d9{x) = 0. 
Jr Jr 
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Let * : M ^ K be such that 

'l ifa;<-A 
otherwise, 



where A > is fixed. Choose = 'i>'{x) (actuaUy one should to use a C°° approximation 

of *'(a;)). Observe that {u^x)x^{uxx) = '^{uxx)x and so Jg^Ux{uxx)x'^' (uxx)d9{x) = 0. Define 
A = {x\uxx < -A}. Thus, 



0= / {ul, + V"{x) + aUxx)Huxx)d9ix) 
(uL + y"{x) + auxx) de{x). 



A 

Since Uxx < —A, and using auxx < — \o? , one can show that, > {\ — + c)9{A), 

where, \V"\ < c. Since A is arbitrary , we get 9{A) = 0. Thus, there exists A > 0, such that, 

Uxx > -A, 61-a.e. 

The solutions of au + ^u"^ + V (x) =0 are semi-concave (this is a standard result, see |BCD97) 
or the survey paper [BG0 8J). so we get that there exists /3 > such that Uxx < /3, and so, for 
some C > 0, \uxx\ < C, 9 almost everywhere. 

4 Holonomic discounted stationary Mather measures 

Motivated by the formal computations in the previous section, we will now establish the existence 
of holonomic discounted stationary Mather measures. In the paper TCG08], these measures were 
called invariant, we did not keep this name here to avoid confusion with invariance with respect 
to Euler- Lagrange equations. 

Given a probability measure v, and a corresponding discounted stationary Mather measure fi 
with trace v, we say that /i is a holonomic discounted stationary Mather measure if 

(/?(0, w)c!/i(u, w) = / ip{0,Lu)di'{uj), 
' X o Jn 

for all ip G Cj^(M" x il). In particular, /i satisfies the undiscounted holonomy constraint. 

Theorem 8. There exists a holonomic discounted stationary Mather measure. 

Proof. Fix Lu Cz fl. Consider a sequence T„ — > cx3 and a sequence Xn{t) of minimizing trajectories 
for the dynamic programing principle (|12p . that is. 



Ua(0,Cj)= f e °'*L{Xn{t),Xn{t),Uj)dt + e "^"Wa(x„(T„),tj). 

Jo 

Because Ua is Lipschitz and 

X„ = -DpH{DxU'^{Xn{t)),Xn{t)) 

the |a;„| is uniformly bounded. 

Define a probability measure ^ by 

1 f^" 

(j){v,uj)dn{v,uj) = hm — / r^„(t)t^)dt, 

for any cj) e C°(M" x fl), where the limit is taken through an appropriate subsequence. This 
sublimit exists and is a probability measure because ^l is compact and \xn\ is uniformly bounded. 
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Let if e Cg. Observe that ^y(a;„(t),a;) = Xn{t)-Dx^{Q,Tx^{t)ijj). So, if (l){v,uj) = v-Dxip{0,u!), 
then 



1 



v,u)d^{v,ijj) = lim — / x„{t) ■ D^Lp{Q,T^^^(^^)Uj)dt 

n^oo In Jq at n^oo i„ 

Since ^"(^ = • Dxip{0,u)) — a(f{0, lo), 

1 r^- 

= hm — / a;„(t) • £'x¥'(0,ra;„(t)w) - a<(C(0,Ta:„(t)w)dt 
n^cx> r„ Jo 

= -a lim ^[ (f{0,T^^(^t)Uj)dt= -a [ (f{0,u!)diy{u)), 

T„->oo In Jo Jn 

where u is given by, 

/ g{w)du{w) = lim ^ [ g{T^^(t)w)dt, 
Jn Tn Jo 

for all continuous function : O — > M. In particular, /gn^o </'(0, w) = Jj^ 97(0, a;)rfi'(a;). 

We must to prove that /i is minimizing. To do so, fix first n and consider a partition {0 = 
to,ti, ...,tN-i = Tn} of [0,T„], where tj+i = ti + h, and h = Tn/N. The restriction of Xn{t) to 
each sub-interval is minimizing, i.e., 

e-'^it-^-) L{Xn{t),X„{t),U>)dt + e~"'^Ua{Xn{ti+l),U!). 

We have, 

i=N-l 

Ma(a;„(ti),w) - e""'*Ma(a;„(ti+i),w) = 

i=0 

i=N-l 

= ^ Ma(a;r»(ti), w) - Ua(a;„(ti+i),a;) + (1 - e""'')ua(a;„(ti+i),a;) = 

i=0 

i=N-l 

ah 

Sending /i ^ we get 

=N-1 

im 



Ma(a;n(0),a;) - Ua(a;„(r),a;) + q;(^— ^^^^ — ) ^ ftua(x„(ii+i), w). 

i=0 



lim V Ua{Xn{ti),Uj)-e~"''Ua{Xn{ti+l),(j) = Ua{Xn{0),Uj)-Ua{Xn{Tn),Oj)+a Ua{Q,T^^(t)^)dt. 

On the other hand, wc have 



lim V / e-"(*-*^)L(x„(t),i„(t),w)di = / L{xn{t),Xn{t),uj)dt. 
i=o Jo 
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Thus, 

a I Ua{0, uj)dv 



rwoo Tn Jq 



lim — / L{xn{t),Xn{t),uo)dt = I L{0,v,uj)dij{v,uj). 



By Corollary [5] we have Ha = a J Ua{0,uj)di'. Thus /i is minimizing. □ 

We should note here that the Theorem does not assert uniqueness. Furthermore the measure 
H may depend on the choice of w € O or in the particular sequence we choose to extract the weak 
limits. For our purposes, however, existence is sufficient. 

Theorem 9. Let fia be a holonomic discounted Mather measure as constructed in theorem HI 
Then /Xq, is invariant under the discounted Euler-Lagrange flow. 

Proof. It suffices to prove that for any bounded function (j>{x,v,uj) G Cj(R" x M" x uj) we have 

W^''\7^^y(t){0,v,u;)dna = 0. 



This follows, from the identity 



(a;„(r„),i„(r„),cj) - (/)(a::„(0),a;„(0),Lj) / —4>{x„{t), Zn{t),^) 



d 





dx dv ' 

dividing by Tn and letting n ^ oo. □ 







5 Graph property, regularity and stationary Mather mea- 
sures 

In this section establish that the discounted Mather measures are supported in a graph of a 
(partially) Lipschitz function. As we are using similar techniques to |EG01| (see also jBGOSJ) we 
will present in this section the main differences and technical points and postpone to Appendix 
Othe detailed proofs. We will the discounted Mather measures to construct a stationary Mather 
measure invariant under the Euler-Lagrange flow. 
We assume that 

L{x + y, V, uj) - L{x, V, uj) <{c + cL)\y\ 

Lemma 10. Let Ua be the viscosity solution of H{Q, DxUa{Q,uj),uj) + aUa{Q,uj) = given by 
Proposition m Then 

lim aua{0, uj) 

a— >0 

does not depend on uj. 

Proof. We know that aua is uniformly bounded, so aua{0, uj) — > ^{uj) pointwise for some function. 
On the other hand, fixed ujq (z fl we know that ?/„(?/, wq) is uniformly Lipschitz in a;, uniformly as 
a — > 0, that is, 

\Ua(xi,UJo) - Ua{x2,UJo)\ < C\xi - X2 | . 
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Thus, if \y\ < R then, 



a— >0 a^O 



that is, lima^o aUa{0,TyUjQ) — hmo^o ^^^(O, cjq) for \y\ < R. 

From Proposition m we know that Wq(0, lu) is Lipschitz in u) with Lipschitz constant K/a^ that 
is, ^ 

|Ua(0,LJi) - Ua{0,UJ2)\ < —d{u!i,UJ2). 

a 

Consider e > and y G R", such that d{TyUJo,uji) < e. Observe that, 

|q;Uq(0,wo) - aUa{Q,uJi)\ < 

< \aUa{0,UJo) - aUa{0,TyLLlo)\ + \aUa{0,TyLLlo) - aUa{0,U!l)\ < 

K 

< \aua{0,ujo) - aUa{0,TyUJo)\ + a—d{TyUJo,uji). 

a 

Sending a ^ 0, and then e — > we get hniQ,^o Q!Wq,(0, wq) = lm'ia^oaua{0,uji). Thus, ^(w) is 
constant. □ 

Lemma 11. Let Ua be the viscosity solution of H{0, DxUa{0,uj),uj) + aua{0,uj) = given by 
Proposition m Then 

hm aUa{0,Lu) = H, 

a— *0 



where 



H = inf I Ldfi, 



and the infimum is taken over aU stationary holonomic measures. 

Proof. Denote by H the hmit as a —> of aua , which is constant by the previous lemma. Let /Zq 
be a holonomic discounted stationary Mather measure. Then, because /Iq is holonomic we have 



H < lim / Ldfia — lim a / Uadfia = H . 

Let be a stationary Mather measure. Then, because ^ is a discounted holonomic measure with 
trace fj, we have 

H — LdjJL > a I Uadfj, H, 

jR^xn JR^xn 

as a ^ 0. This shows that H = H. □ 

Lemma 12. Let fia be a sequence of discounted stationary Mather measures with trace i>a- 
Suppose that /i^ — > /i when a 0, then /i is a stationary Mather measure. 

Proof. First we must to prove that /i is a holonomic probability measure. In fact, for any ip G C^, 

V ■ Dxip{0,uj)diia — / (p{0,uj)dfj,a — a / ip{0,uj)dva ^ 0, 
"xo JR^xn Jn 

when a ^ 0. 

Using Corollary we get 

/ Ldfj, = lim / Ldfia = lim Ha — lim / aUa{Q,uj)di'a{^) = H. 

Thus is a Mather measure. □ 
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Theorem 13. Let /ic, be a discounted Mather measure with trace i^a (or if a = a stationary 
Mather measure). Then /i^ is supported in a graph, that is, there exists a measurable function 
Va-.n^R" such that, 

supp^Q = {{v,uj) e R" X n\v = Va{l^)}. 

Proof. As in |BG08| . for instance, we just observe that the result follows from the fact that the 
Lagrangian is strictly convex in v, whereas the discounted holonomy constraint is linear in v. □ 

Since the holonomic discounted measures are also holonomic, the same techniques in |EG01| 
(see also |BG08p can be adapted to establish the following regularity result: 

Theorem 14. Let Ha be a holonomic discounted Mather measure. If Ua is a viscosity solution 
of Uni), then for each y £ M, 

\DxUa{y, uj) - Dj;Ua{0, uj)\ < C\y\, 

9 almost everywhere and uniformly in a. 

The proof of this theorem since it follows (almost) exactly the same steps as in |EG01j (see 
also [BG08 ) and is presented for completeness in appendix [Cl The only difference is the term au 
in the Hamilton- Jacobi equation, which can be controlled, as discussed in section [3l because we 
are using holonomic discounted measures. As a corollaries to the previous Theorem we have 

Corollary 15. Let be a holonomic discounted Mather measure. Then, there exists a function 
: n — > K", such that supp/ic — {{v,uj) E R" x il\v = Va{u!)}- Furthermore, Va is partially 
Lipschitz in the following sense: 

\Va,{TyLu) - Vc.iu)\ < C\yl 
for all Lu in the support of fia, and C is uniformly bounded as a ^ 0. 

Finally, our last result concerns the existence of stationary Mather measures invariant under 
the Euler-Lagrange flow. 

Theorem 16. There exists a stationary Mather measure /i which is invariant under the Euler- 
Lagrange flow. Furthermore /i is supported on a graph. 

Proof. Let fi^ be holonomic discounted Mather measures as constructed in theorem [S] Consider 
a weak limit /z. By lemma [T^ /i is a stationary Mather measure. Because for any (j){x, v) we have 



I 



W"^" V,,,(/.(0, V, uj)dfi^ = a / {Dt^L)-'D^LD^(t){Q, v, w)dMa. 

xSl JE'>xO 



we conclude that 

/ VF^°V,,,0(O,w,w)dM = O. 
jR"xn 

The graph property of stationary Mather measures follows from theorem 1131 □ 

A Proof of Theorem [1] 

In this appendix we present the proof of Theorem [2 as well as some background material. 
Let 7 be as in Let A4 be the set of weighted Radon measures on 17 x M", i.e., 

M — {signed measures on M" x with / ld\fi\ < oo}. 

JR" xSl 
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Note that M is the dual of the set CiJ(R" x Q). 
We introduce the foUowing sets 

A^i = {/iGA^|/iisa positive probabiUty measure}, 

and 

M2 = {f^eM\ I V- D^ip{0,uj)dn{v,uj) = 0, for all ip £ C](R" x n)}. 
Using this notation the Mather problem can be reformulated as 

min / L{0,v,LLj)dii{v,uj). 



MinM2 



Consider the following subset of functions : M" x ^ M, 

C = c\{(f> I (l){v,uj) = V ■ D^<f{0,uj), for some Lp G Cl{W x Q)}. 

Observe that C is a closed convex set. 
For (/> G C° let 

h{(l3) = sup uj) - L(0, V, Lu)). (15) 

Since h is the supremum of linear functions, it is a convex function on C^. As we will see bellow 
in Lemma [T51 ft, is a continuous function. 



For </)G C^, let 



nw = i ^, . (16) 

— cx) otherwise. 



As C is a closed convex set we have that g is a concave and upper semicontinuous function. 
Therefore its Legendre-Fenchel transform is given by 



g*{^i) = inf - / (j3d^i ~ g{(l)) . (17) 
Since /i is a convex and lower semicontinuous function, its Legendre-Fenchel transform is given by 



h*{^x)= sup -/ (t)diJL ~ h{(t)) ] . (18) 



Proposition 17. Let 5 and ft defined as in (|T5|) and (fT6|) . Then 



and 



otherwise, 



if /i G A^2 
— cx) otherwise. 



Proof. First we assume that G A^i. As ft is a convex function, its Legendre transform is given 
by dH]). Using the definition of ft, equation (fTKjl . we get 



ft*(/i) = sup (— / — sup (— 0(u, w) — L(0, u, w)) 

0eC«(R"xn) V JR"xn R"xf2 

Consider the family of compact subsets of R" x Vt given by 

Kn = {(w,u;) G R" X f]| |w| < n}. 
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and let i]n : R" x i7 ^ R be a continuous function such that < 77„ < 1, = 1 in Kn, and 
suppr^n C Kn+i- Then define 

Observe that the sequence L„ is increasing and pointwise convergent to L{0,v,lu). 

Is easy to see that L„ G C°(R" x il). Furthermore, for fixed n, one can write any function 
e C°(R" X fl) as (j) = —Ln + "0 where V' & C°(R" x Q). From this observation we get 



= sup - / {-Ln + ^|J)d^l- swp {-{-Ln + Tp) - L) 



sup I / Lndii — / ^pdiJ — sup [Ln — L — ip) 

ipec°(R"xn) V"'R"xo Jr-xh r^xO y 

/ Lnrf/^ + sup / I ~1p — sup {Ln — L — -0) I C^M- (19) 

■/r"xo V6C'2(R"xn) jR"xn V R"xn / 



)/'eci;(R"xn) JR"xn 

If we take = in (fT^ we have 



/ i„d^ + / I - sup (L„ - i) I d/i > / Lndfi. 

JRixf2 JR"xSl V R"xf2 / jR'>xn 



Thus using the monotone convergence theorem we get 



h*{fi) > [ Ldfi. 
jR^xn 



'xSl 

In order to get the other inequality we can rewrite (fT9|) as follows 



= / Lndfi + sup / I {S — ^ — S) — sup {S — tJj) \ d/i 
JR^xn i/;eci; JR"xn V r-xo / 

— / Lndfi + sup / i {S — ip) — sup (S" — ^) j d/x — / Sdfi, 

where S = Ln~ L. Since /i e A^i, we have /r^xo (("^ ~ 0) ~ supRnxfil"^ ~ "0)) '^M 0- Therefore 

< / Lndfi — / Sdfi = / Ld/i. 

^R^xO JR"xa JR^xn 

If /X ^ A^i, we have two possibilities. First, if ^ then we can find a positive function 
ip G C°(R" x ^1) such that / ^pdfi < 0. Define i^n = ni^ G C°(R" x ^2), then 

h*{fJ.)> I - / V'nC^M - sup {-tpn - 



-ipdn + inf {ip H — L) ) +oo. 



'xSl 



'xn 



when n ~^ oo. 

On the other hand, if /i > but J dfi ^ 1 we take = A; e R, then 



/ kdfi — sup (— fc — L) j 

jR"xO R"xO / 



fc 1 - / + inf L 



'xn 



'xO 
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when k — > ±cxd, because L > 0. 

Now we compute the Legendre transform of 5. As g is concave we compute its Legendre- 
Fenchel transform using (fT7|) . First we suppose /i e Ai2- In this case we have two possibihties, if 
(j) eC then 



else if, ^ C then 



(pdn - g{4>) = 0, 

"Xf2 



— — — I (j)dfi — (—00) = +00 



thus = 0. 

Otherwise, ii fJ, ^ A^2 there exists 4){v, us) — v- Dx(p{0, w) € C such that /u„xf2 '^'^/^ ^- Define 
4>\ — \v ■ Dxip{0, w) e C then 

9*i^J■)<{- 4>xdfi- g{(l3x)] = (j)d^ -^-CG 

when A — > +00. □ 
Remark 3. Observe that 



In fact. 



Lemma 18. The function 



is continuous. 



' j L{(),v,uj)dii{v,uj) ~ if ^leMlC^M2 

J L{0,v,uj)dfi(v,Lu) — {~oo) if fi e Mi a,nd fi ^ M2 

+00 ~ (0) if fi 1^ Ml and /i e A^2 

+00 — (— 00) if fi ^ MiSLiidfi ^ M2- 



h{(f)) = sup cj) — L(0, ti, Lj)) 

R"xO 



Proof. Let (/)o be an arbitrary, but fixed, function in C^(R" x il). Suppose (/)„ 0o, that is \\(j>n — 
4'o\\i ^ 0. Let B^{(j)f)) = |0 G C°(M" X Vl) \ ||0„ — 0||-y < e} be the baU of radius e centered in ^q. 

Take e 5e(<^o)- Since, hm|„|^oo = \mi\^,\^^ ^'^^(v)^ = +°° ^nd \im\^\^^ '"^"^^"j"^ 
uniformly on w € fi, given S,M > 0, there exists i? > such that 

l%^l<<5 if \v\ >R 
i^>M if \v\ >R. 

Then, for > i?, 

N \ f -(t){v,uj) + (j)o{v,uj) -(j)o{v,uj) L{0,v,uj) \ 

-(^{v,uj) ~ L{0,v,uj) = < — + — — < 

I 71^') 7W 7(w) J 

< ^{£ + 3- M)\v\ ^ -00 
\v\ 

when \v\ +00. 
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As 

\\<P,i - (Poh = sup —■ > 

K"xf2 

we have that, for n big enough, wc can choose R in such way that 

h{(f)) = sup (— w) — L(0, w, w)), 

{\v\<R}xn 

and 

h{(f>n) ^ sup (-(^„(w,u;) - L(0, 

Since the convergence —<pn — L{0,v,lu) — > — L(0,v,u!) is uniform on the compact {\v\ < 
R} X n, we have 

hm h{(f>n) — hm sup (— (^„(u, oj) — X(0, w)) = 

„^oo "-'°°{|t,|<fl,}xJ2 

= sup hm (-0„(w,cj) - L(0,w, w)) = /i((/)o). 

{|t;|</?}xn"-°° 

Thus the lemma is proved. □ 

The last ingredient of the duality is the Legendre-Fenchel-Rockafellar Theorem, see for instance 
[Vil03) . 

Theorem 19. (Legendre-Fenchel-Rockafellar) Let be a locally convex Hausdorff topological 
vectorial space over M with dual E*. Suppose that h : E ~* (— cx3,+cx3] is convex and lower 
semicontinuous and g : E* —^ [—00, +00) is concave and upper semicontinuous. Then 

min(/i* — g*) = sup((7 — h), 

E' ' E 

provided that /i or g is continuous at some point where both functions are finite. It is part of the 
theorem that the left hand side is a minimum. 

Lemma 20. Define the functional, S^cf)) = g{4)) — h{(j>). Then S is uniformly continuous in the 
interior of C. 

Proof. In fact, given e > 0, if — (/'2II7 < e, that is, —ej{v) < (j)i{v,w) — 4'2{v,w) < sj{v), for 
all {v, w), then 

\SiM~Siq^2)\<{infjiv)}e. 

In particular 

sup g{(t)) - h{(j)) = sup g{(t))-h{(t)). 

'■peel 

□ 



B Proof of Proposition [4] 

Proof oj Proposition^ We must to prove that the function Ua is stationary. Since L > Q, Ua 
is well defined as an infimum. On the other hand the stationarity is an easy consequence of 
the correspondence between the set of all globally Lipschitz trajectories with initial condition 
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a;(0) = X and the set of all globally Lipschitz trajectories with initial condition y{0) — 0, given by, 
{x{t)} {y{t) = x{t) - x}. In fact, 



at 1 



Wa(0,T^w)= inf / e "*L{y{t),y{t),T^Lu)dt 
a(o)=o Jo 

r+oo 

= inf / e~°'*'L{{x{t) — x) + x,x{t),Lu)dt — Ua{x,Lu). 

In order to prove that Ua is a viscosity solution in w, let 1^9 : R" x J7 K be a stationary 
function such that Uq,(0,w) — ip{0,Lu) has a local minimum (resp. maximum) in lu^ e il and 

Ua{0,UJ^) - ip{0,UJ^) = 0. 

Consider a trajectory satisfying x{0) — such that x{t) is a finite time minimizing, globally 
Lipschitz trajectory, for the dynamic programing principle (HH), that is, 

Wa(0,Wy)= / e-"'L{x{t),x{t),uj^)dt + e-"^Uo,{x{T),uj^), (20) 
Jo 

for T small enough. 

Suppose that Tl(ip,uj^) < 0, by continuity there is a neighborhood B of uj^ in and ^ > 
such that H{(p,uj) < —S for all uj £ B. Since Ti.{(p,uj) — H{0, Dx(p{0,Lo),uj) + a(p{0,uj) we have 
—vDxfiOTCu) — L{0,v,uj) + aip{0,uj) < ~S, for all w G i? and v £ M". If we choose v = x{t) and 
to = Tj-j-j-jWip then 

x{t)D^ip{0,T^(t)^^) + L{x{t),x{t),uj^) ~ a(p{x{t),uj^) > S, 

for < t < T. 

Integrating this expression and using ^<p(x(i), w) = x(t)Dxip{0, T,j.(^t^uj) we get, 

'p{0,Tx{T)i^ip) - ^iO,^ip) + / L{x{t),x{t),uj^)dt ~ a ip{x{t),Lu^)dt > ST. 

Jo Jo 

Since Ua{0,u!) > (^(0,w) in B and UQ,(0,a;i^) = 1^9(0, W;^), we have 

Ua{0,Tx{T)^^) - Ua{0,UJ^) + L{x{t) , x{t) , uj^)dt - a ip{x{t),uj^)dt > ST. 

Jo Jo 

Using ([201) in the last inequality we get, 

(1 - e-"^)u„(0, T^(T)^^) + [ (1 - e-°'')L{x{t), x{t),uo^)dt -a f ^{x{t), uj^)dt > ST. 

Jo Jo 

Writing 

T i 

Ua{0, T^iT)UJ<p) + ^ g-aT-) y (1 - e""*)L(x(t) , i (i) , w<^)di- 
T I F T 

and usmg hin ^ _ ^_^t 

"a(0,Wip) - (^9(0,0;^) > — 

contradicting Ua{Q,ijJtp) = ^p{Q,uJ^p). 

The proof for the maximum case is analogous and so the theorem is proved. □ 
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C Proof of Theorem [14 



In this last appendix we give a proof of Theorem 1141 Before that we need to estabhsh some 
additional results. We note here that we will be using the techniques in [EGOl] (see also |BG08| ) 
adapted to the stationary setting. 

Remark 4. Let Ua be a viscosity solution in uj of (fTU|) . then, because it is also a viscosity solution 
in X (Proposition [2]) and it is Lipschitz, DxUa{0,TyUj) is defined Lebesgue almost everywhere and 

H{y,D^Ua{y,uj),uj) + aua{y,uj) = 0, 

for Lebesgue almost everywhere y £ R". 

For any probability measure we can define a new measure of probability p, in R" x D, given 

by, 

/ ^{Pi^)dfL{j),uj)^ I ■>p{-DyL{0,v,uj),u)dfj.{v,u). 
In this case, the integral holonomy constraint can be rewritten as 

DpH{0, p, Lu) ■ DM0, ^)dfi{P, OJ) = 0, 

Vv?e ci(R" X n). 

Theorem 21. Let fia be a holonomic discounted stationary Mather measure. Denote the projec- 
tion in the coordinate uj of iia by 9a, that is 

(p{uj)d9a — / tp{u!)dlJa- 



If Ua is a viscosity solution of (|10p . then DxUa{0,(-u) exists 0Q.-a.e, and /ia-a.e, p = —DxUa{0,Lu). 

Proof. By the strict uniform continuity of H there exists 7 > such that for any p,q,y (z R" and 
w g 51 we have 

H{0,p, Tyu) > H{0, q, TyLo) + DpH{0, q, Tyu){p - g) + - gp. 

Let u"^ = u * 77, by Remark^ for almost every uo and y, let p — DxUa{0, Tyu) and q = D^u'^aiO, 
Then 

H{0,D^Ua{0,TyUj),TyUj) > H{0,D^u''JO,Uj),TyUj) 

+ DpH{0,D,U%{0,Uj),TyUj)iD,Ua{0,TyUj) - D.,U%{0,UJ)) + ^\D,Ua{0,TyUj) - D.,U%{0,UJ)\^. 

Multiplying by 77"^ (y) and integrating we get 

H{0,D,uU0,Lu),Tyu)7j%y)dy+ f ^\D^UaiO,TyU;) - D,uliO,iu)\^ii%y)dy 



< / H{0, DxUa{0,TyUj),TyUj)ri''{y)dy 

DpH{0, D,u%{0, uj), TyLj) [D,u%{0, uj) - D,Ua{0, TyCj)] Tj'{y)dy. 
Remark HI implies that, H{y, D^Uaiy,'-^),'^) — ~OLUa{y,to) almost everywhere y. Thus 

H{Q,D,ul{Q,Tyu),Tyu)Ti^{y)dy + (3'{w) < -a<(0, c^) + o^(£) (21) 
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where 



On the other hand, the convexity of H, imphes that, 
l;\Dxul^{Q,uj) - p\^djla{p,uj) 
< f [HiO,Dxul^{0,iu),Lu) - H{0,p,Lu) - DpH{0,p,io)iDxuliO,u;) -p)]dila.ip,uj) 
H{0, D^ul^^O, uj), uj)dp.a{p, uj) 

[H{0,p, uj) + DpH{0, p, uj)D,u% (0, uj) - DpH{0,p, iu)p] dfi^ {p, oj) 

'xf2 

H{0, D,ul{0, u), u) + L(0, ^DpH{0,p, Lo),Lo)dflo,{p, to) 

= / H{0,D,u%{0,uj),uj)dtlc,{p,uj) + H^. (22) 
jR"xn 

Integrating (PT|) with respect to fl and adding (|22p . we get 

^\D,ul{0,Lu)-p\^dfia{p,uj)+ [ l3%Lu)de{uj) <o{e). 
So, 9a almost everywhere we have DxUa{0,uj) — hm Dxu'^a{0,uj), in particular p = DxUa(0,uj) in 

e— »0 

the support of //. □ 

Theorem 22. Let be a holonomic Mather measure for the discounted stationary Mather 
problem. If Ua is a viscosity solution of PT7| . then for each ft, e M, 

almost everywhere. 

Proof. If ft 7^ then we define, 

Uaix, Uj) = Ua{x + ft, Uj) and itaix, uj) — Ua{x — ft, Ll)), 

and u^(a;,a;) and {{^(Xjw), the corresponding smoothings (see Remark [1]). 
Remember that 

H{h, Dxu%{0, uj),Lo) + au'jO, u) < ce, 

and 

H{~h, D^u^'aiO, Lu),Lu)+ au^(0, uj) < ce. 

Thus, 

iJ(0,D,w^(0, uj),uj)^ 2H{0, D,ua{0, uj), uj) + H{0, D.u'jO, uj), uj) 

=H{0, D,u%{0, uj), uj) - H{h, -D,<(0, io),uj)+ H{h, D,ul{{), uj), lu) + aul{Q, uj) 
- (0, w) + 2aua (0, w) - au^ (0, w) 

+ a?i^(0,w) + iJ(-ft,i?xU^(0,w),cc>)-i/(-ft,i^,w^(0,w),w)+iJ(0,^,u^(0,tj),tj) 
<2ce- a{u%{0,ui) - 2u„(0, w) + w^(0, w)) 

-{DxH{0,Dxul{0,Lo),uj)-DxH{0,D,uliO,u),u))h + Oi\h\''). (23) 
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On the other hand the convexity of H imphes that 

i7(0, Z?.<(0, oj),oj) >H{0, D,u^{0, uj), uj) 

+ DpH{0, D,u„(0, uj), uj){D^ii'^{0, w) - -D,u„(0, tu)) 

+ ^\D,uUO,uj)-D,u^{0,oj)\^, 



and 



77(0, D,u%{0, iu),iu) >H{0, D^u^{0, u), 

+ DpH{0, D^u^{Q, Lo), c^)(I?,u^(0, oj) - D.,u^{Q, lo)) 

Adding these two formulas we obtain the foUowing inequahty: 

I {\D,ul{Q,uj) - D,w„(0,c^)p + \D^ul{Q,u) - D,u^{Q,uo)\^) 

+ DpH{Q, D,uo,{0, uj), uj){D^ii'^{0, uj) - 2D,u^{0, uj) + D^u^iO, uj)) 

< {H{0, D^uUO, uj),uj)~ 2H{0, D,Ua{0, uj),uj)+ H{0, D,u%{0, uj), uj)) . 

By (^51) we have, 

I {\D,ul{0,uj) - D,ua{0,uj)\^ + \D,ul{0,uj) - D^u^{^,uj)f) 

+ D.pH{^,D,Uo,{^,uj),uj){D^ul{^,uj) - 2D,u^{^,uj) + D^u'J^^,uj)) 
< 2c£- a(M^(0,w) - 2ua{0,uj) +ul^{0,uj)) 

- iD,H{0, D^iiliO, uj), uj) - D,H{0, D.uKO, u), u)) h + 0{\h\^). 
Or equivalently, 

^ {\D,ul{0,uj) - D,u^{0,uj)\^ + \D,u%{0,uj) - D,u^{0,uj)\^) 

+ DpH{0, D,u^{0, uj), uj){D^ul{0, uj) - 2D,u^{Q, uj) + D,ul{0, uj)) 
+ a {ul (0, uj) - 2u^ (0, uj) + ul{Q,uj)) 

< 2ce - {D,H{0, D^iiliO, uj), uj) - D,H{0, D^u^iO, iu),iu)) [h) + 0{\h\''). 
Define, j3'^{x,uj) = u%{x,uj) — 2ua{x,uj) + u'^{x , uj) , so can be rewritten as 



^ {\D,ul{0,uj) - D.,uc,{0,uj)\^ + \D,ul{0,uj) - D.,u^{0,uj)\' 
+ DpH{Q, D,uo.{0, uj),uj)D^(3'{Q, uj) + a/3^(0, uj) 



< 2ce - (D^HiO, D^niiO, uj), uj) ~ D,H{0, D.uKO, uj),uj)) [h) + 0{\h\'') 

Applying the inequality 

II {D,H{0, D^uliQ, uj), uj) - D,H{Q, D,ul{Q, uj), uj)) h\\ 
<\\D^pxH{Q,D,ul{Q,uj),uj)\\ ■ \D,ul{Q,uj) - D,ul^{Q,uj)\ ■ \h\ 

<^\D^ul^{Q,uj) - D^ul{Q,uj)\'' + -\D^uUQ,uj) - D^uU^,uj)\'' ■ \h\^ 
4 7 

<\ {\Dxul{Q,uj) ~ i?,u„(0,w)|2 + \D^ul{^,uj) - i?,u„(0,c^)|2) 

+ -\D.,u^^(S),uj)^D,if^(S),uj)\^ ■\h\\ 
1 
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to ([25]) we get, 



+ DpH{0,D^Uc,iO,uj),uj)D^(3^0,Lu) + aP^iO^Lu) < 2C{e + \h\^). (26) 



Consider a function * : M ^ M, such that $(s) = *'(s) > 0. We can muhiply ([Ml) by $ (^^Tqr^) 
and integrate with respect to 0, 



+ I DpH{0,D,u^{0,u;),u;)D,P^O,uj)<^> { ^^^^^ ) dO 



We have 



\h\^ I^DpH{0,D,u^{0,u;),u;)-^D,P'{0,u;)^ (^^^^^ ) 
DpH{0,D,u^{0,u),u)D.,^ (^^^ {0,uj)de = 0. 
Thus, (|27p can be restated as. 



f \D,uliO,iu)-D,uUO,iu)\^^ 
Jn 



dO 



Define, A\ — {a;| ^ i^^a"^ < —A}, and consider the function ^ defined by 



h|2 



s if s < —A, 
1 otherwise. 



Fix a positive constant 7 such that the functions — Ua{x,uj) — ■^\x\'^ and u^^{x,uj) — 

u^{x,u!) — ■^\x\'^ are concave. Observe that a point uj is in A\ only if 

u^(/i,t^)-2u„(0,t^) + u^(-/i,a;) < -(A + 7)1^1^- 
Define F'^{t) ~ u^(t|^,a;). Since F"^ is concave and (F^)" < we have 

c^) - 2u%{0, Lu) + uli-h, uj) > [D^uliK ^) - D,ul{-h, Lo))h. 
Subtracting this incquahties we get, 

(A + 7)|/ip < 2|i2^(0, Lu) - w„(0, + \D,ul{K uj) - D,uU-h, \h\. 

Since Uq is stationary and uniformly Lipschitz continuous wc have \u'^^{0,lu) — Ua(0,LL!)\ < Ce. 
thus we can choose s in such way that 

iD^u'^ih.io) - D,uU-h,u;)\ > i^+j)\h\ 
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and 

\D,ul^{h,u;) - D,ul^i~h,u;)\ > i^+^)\h\. 
Using this estimates in (6) we get 

+ l)'\h\^0iAx) - aX\h\^0{A^) < 2C{e + 

Observe that, if 6{A\) > then the left hand side of this inequahty converges to +00 when 
A —)■ +00, so there exists a value Ao such that 0{Ax) — 0, that is, — Ao|^P < u'^ai^i ^) ~ '^Ua{x, uj) + 
'u^(a;,a;), almost everywhere. The upper bound comes from the semiconcavity oiua. Thus there 
exists C > such that \ua{h,uj) — 2ua{0,uj) + Ua{~h,uj)\ < C\h\'^, 6 almost everywhere, which 
completes the proof of the theorem. □ 

Proof of Theorem\14\ Let 9 be the projection of /x^. By Theorem!^ DxUa{0,uj) exists d-a.e. On 
the other hand, fixed uj G supp^, DxUa{y,uj) exists Lebesgue almost everywhere. 
We claim that 

\Ua{y,Uj) - Ua{0,Uj) + Ua{~y,uj)\ < C\y\'^ . 

This claim is a consequence of Theorem [22l by choosing h = y and of the semi-concavity of Ua ■ 
In fact, we have 

- C\h\^ < u^{y,to) - 2uo.iO,L^)+Uo.{~y,Lo) < C\h\^, (29) 
Uaiy, uj) - Mq(0, u) - D^UaiO, io)y < C|?/p, (30) 

and 

Ua{-y, uj) - Ua{0, uj) + D^UaiO, Lu)y < C|y|^. (31) 

The claim is obtained from (PU]) and from the difference between ([^H) and (pij) . 
Let z G M be a point such that \z\ < 2\y\. The semi-concavity of Ua implies that, 

Ua{z,uj) < Ua{y,uj) + D^Ua{y,uj){z -y) + C\z - yp. (32) 

Using, Ua{z,u) = Ua{Q,uj) + DxUa{Q,ui)z + o{\z\'^) and Ua{y,uo) = Ua{0,oj) + DxUa{0,ijj)y + o{\y\'^) 
in (4) we get 

(Z?,u„(0, Lo) - D^uo^iy, ~y)< C\y\^- (33) 
If we take z = y+\y\ \Dlll^(a%]~Dlll\lZ]\ then we obtain \D^Ua{y,uj) - D^Ua{0,uj)\ <C\y\. □ 
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